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We present a new strategy for measuring the electron neutrino mass (m Ve ) by future detection 
of a Galactic supernova in large underground detectors such as the Super-Kamiokande (SK). This 
method is nearly model-independent and one can get a mass constraint in a straightforward way from 
experimental data without specifying any model parameters for profiles of supernova neutrinos. We 
have tested this method using virtual data generated from a numerical model of supernova neutrino 
emission by realistic Monte-Carlo simulations of the SK detection. It is shown that this method is 
sensitive to m„ c of ~ 3 eV for a Galactic supernova, and this range is as low as the prediction of 
the cold+hot dark matter scenario with a nearly degenerate mass hierarchy of neutrinos, which is 
consistent with the current observations of solar and atmospheric neutrino anomalies and density 
fluctuations in the universe. 

PACS number(s): 95.55.Vj, 14.60.Pq, 95.35.+d, 97.60.Bw 



I. INTRODUCTION 

It is well known that detection of a neutrino burst from 
a collapse-driven supernova by large underground detec- 
tors gives us some constraints on neutrino masses, due 
to the delay of arrival times depending on the neutrino 
energy as At = 5.15(Z?/10kpc)(r7i l/ /leV) 2 (e t/ /10MeV)- 2 
msec, where D is a distance to the supernova, m„ the 
neutrino mass, and e„ the neutrino energy |IJ. Some up- 
per bounds on the electron neutrino mass (?n„ c ) , ranging 
in 10-20 eV, have already been derived by a number of 
papers using the historical data of SN1987A On the 
other hand, given the status of tritium /3-decay experi- 
ments, experimental upper limit on m„ t are also consid- 
ered to be ~ 10-15 eV j|. Therefore the next Galac- 
tic supernova expected in the near future and its detec- 
tion by the currently emerging international network of 
second-generation neutrino detectors, such as the Super- 
Kamiokande (SK) Q| or SNO ||, would give an impor- 
tant opportunity of setting a more stringent constraint 
on m„ e than the current astrophysical or experimental 
limits. Especially, the normal water Cerenkov detectors 
are the most sensitive to electron antineutrinos (y e 's) 
and an enormous number of v e p — > e + n events expected 
in the SK (~ 5000-10000 events) for a supernova at the 
Galactic center (D = 10 kpc) would give much better 
statistics than that of SN1987A. 

However, the electron neutrino mass measurement by 
supernova neutrinos generally suffers significant uncer- 
tainties related to the original profiles of supernova neu- 
trino emission, i.e., neutrino luminosity curve and energy 
spectrum. The analyses on SN1987A data were based on 
the luminosity decay during the cooling phase of hot neu- 
tron stars, during which the majority of D^s is emitted. 
Because the decay time scale in this phase is O(10) sec- 



onds, it is difficult to probe the arrival time delay shorter 
than this scale. This is why we could not probe the mass 
scale smaller than - 10 eV for SN1987A (D = 50 kpc). 
Therefore it is clear that we have to devise a different 
strategy which maximally utilizes the much larger num- 
ber of expected events in the SK. Although the cooling 
phase is ~ 10 seconds long, the time scale of initial rise of 
neutrino luminosity is much shorter: a recent numerical 
simulation of gravitational collapse and neutrino emis- 
sion H shows that this time scale for zVs is 1-10 msec. 
The prediction of this time scale by the current theory of 
collapse-driven supernovae is robust because this is de- 
termined by the time scale for the shock wave generated 
by the core bounce to cross the neutrino sphere j?J . This 
suggests that we can probe the neutrino mass of ~ 1 
eV, at least in principle, by using a sufficient number of 
events around the initial steep rise of neutrino luminosity. 
In the following we propose a new strategy to set a con- 
straint on m„ e from the initial rise of v e p events assuming 
a detection by the SK. We then test this method by vir- 
tual data of neutrino events detected by the SK, which 
are produced by realistic Monte-Carlo simulations (MCs) 
with a numerical model of supernova neutrino emission 
||. We find that m„ c of ~ 3 eV can be probed by a 
future Galactic supernova. 

II. GETTING CONSTRAINTS ON THE 
ELECTRON NEUTRINO MASS 

Strong v e emission suddenly breaks out when the shock 
wave passes the neutrino sphere with a time scale of 1- 
10 msec, and after this breakout the time variability of 
neutrino luminosity or energy spectrum is on a scale of 
~ 1 second 0-0] . The signature of a finite neutrino mass 
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which we try to detect is the earlier arrival of high en- 
ergy neutrinos in the breakout. Since events from the 
v e p — > e + n reaction are dominant, we treat all events 
as this reaction, and the validity of this approximation 
will be checked later. Suppose that we get a sequence of 
arrival time and detection energy of positrons as (£1, £i), 
(^,£2), • • •) (ijVjEjv); where N is the observed number 
of events, and order of events is defined as tk < tk+i- 
Consider a transformation of detection time of events 
(tk — > t' k ) for a given value of m v< ,, which subtracts the 
arrival time delay due to the assumed neutrino mass, as 



Then we get a constraint on m Ve with n sigma confidence 
level, for a given value of N cut as 



t' k = tk 



Dm 2 



2c(e k 



(1) 



where A np is the neutron- proton mass difference (= 1.3 
MeV). Let e' k be the sequence of detection energy in 
increasing order of t' k . Since the neutrino spectrum is 
roughly constant after the breakout on a time scale of ~ 
1 second, it is expected that the distribution of e' k is ran- 
dom without any correlation to t' k , if the assumed value 
of rn v is correct. Here we define a measure of correlation 



between t' k and e' k in the first N cut events as follows: 



S\m Va ) = Y. 



k=2 
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where Nk(e k ) is the number of events detected earlier 
than t' k with energy greater than e' k , and /(e) is the frac- 
tion of expected events with energy greater than e. We 
can calculate /(e) from v e spectrum, the cross section 
of 9 e p reaction, and detection efficiency of the SK. If we 
use the Fermi-Dirac (FD) distribution with zero chemical 
potential as the spectrum of neutrinos, we can straight- 
forwardly calculate S 2 as a function of m Va , FD temper- 
ature T Pe , and iV cu t from a given experimental data set 
of (tk,£fc). It is known that the real energy spectrum 
of supernova neutrinos is slightly different from the pure 
black body radiation because of energy-dependent opac- 
ity of neutrinos. However we note that we are paying 
attention to how random the e' k distribution is, and a 
detailed shape of the spectrum is not important in our 
analysis. This will be checked later. Now let us consider 
the physical meaning of S (m l , c , T Pe , N cut ). If we approx- 
imately regard the Poisson distribution as the Gaussian 
distribution, the distribution of S 2 is the x 2 statistics 
with iV cu t — 1 (~ N cu t) degrees of freedom, and hence 
it is expected that S 2 obeys the x 2 distribution if the 
assumed m„ c is correct. On the other hand, if the as- 
sumed is significantly different from the true value, 
only high- or low-energy neutrinos will arrive earlier and 
S 2 will become larger than the expectation from the \ 2 
distribution. Similarly, incorrect values of Tp e will lead to 
unexpectedly large S 2 . Therefore S 2 is expected to take 
the minimum at the correct values of m„ c and T Pe . Here- 
after we always take a value of Tp e which minimizes S 2 . 



JVcut) + ny/2N cut , (3) 

with the best-fit value of m„(N cui ) which minimizes S 2 , 
where 27V cut is the variance of the % 2 distribution with 
TVcut degrees of freedom. We stress that this strategy 
does not require any specification of model parameters, 
and constraints on m Ve can be calculated in a straight- 
forward way from experimental data. It should also be 
noted that we have implicitly assumed in the above argu- 
ment that the distance to a supernova is known, although 
in a future detection it may be unknown. In such case 
the above strategy is still applicable but we can only get 
constraint on Dm 2 . 



III. TEST FOR THE STRATEGY BY 
MONTE-CARLO SIMULATIONS 

In the following of this letter, we give a test about the 
reliability of this strategy, by using neutrino emission pro- 
files of a one-dimensional numerical model of supernova 
explosion which does not assume any particular energy 
distribution of neutrinos ||. This is a model of SN1987A 
with a main-sequence mass of ~ 20M Q . We have made 
virtual data sets of by Monte-Carlo simulations 

of the SK detector supposing a supernova at a distance 
of D = 10 kpc. The MC simulation is described in Ref. 
B which takes account of the SK detection efficiency, en- 
ergy resolution, and reaction modes of D e p absorption, ve 
scatterings, and charged-current v e (v e ) absorptions into 
oxygen. Therefore we can check the validity of the ap- 
proximations in the proposed method, i.e., assuming the 
FD distribution and regarding all events as v e p events. 
We have made four MC realizations with simulated neu- 
trino masses (m^ ) of 0, 3, 5, and 7 eV, and applied the 
method to these data varying the value of AT cut . We have 
used events with detection energy greater than 10 MeV, 
to avoid the background noise of (v, v'pj) and (y, 1/717) 
reactions on 16 O ||. The obtained best-fit m^* and 2 
sigma (95% C.L.) lower- and upper-limits (rn l v and m™ ) 
are shown in Fig. 1 as functions of Af cut . If the strategy 
correctly detects the signature of a finite neutrino mass, 
the best- fit mj' should not vary with 7V cut . In other 
words, the constancy of mj* against N cut gives an impor- 
tant consistency check of this analysis. The figure shows 
that, for m^ LC = 3, 5, and 7 eV, is almost constant 
at the simulated values in A^ cut ~ 100-300, suggesting 
that this method correctly detects a finite neutrino mass 
if we use neutrinos of the first 200-300 events, i.e., during 
60-80 msec after the core bounce. One can see a slight 
systematic decrease of mj' in iV cu t > 300, and this is 
an effect of spectral hardening of v e 's, which begins from 
~ 100 msec after the core bounce as a signature of the 
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delayed explosion mechanism It should be noted 

that, in the prompt explosion scenario, stellar envelope 
is expelled in a very short time scale of 1-10 msec which 
is the same as that of the shock breakout, and neutrino 
spectrum is rather constant after the breakout. Therefore 
the proposed strategy would work even better in prompt 
explosions. The optimal value of N cnt depends on the 
distance and profiles of neutrino emission, and it should 
be determined by the constancy in the m^*-iV cu t diagram 
produced from real data detected in the future. 

We next estimate the sensitivity of the proposed strat- 
egy by statistical average of many MC realizations. One 
hundred MC realizations are generated from the super- 
nova model with D = 10 kpc for each of three values 
of m^ c = 0, 3, and 5 eV, and the proposed method is 
applied to these data with iV cu t = 200. The average of 
ra^raj, , and m" for the 100 MCs are shown in Table 1 
with 1 a statistical fluctuations. In order to check the va- 
lidity of estimated confidence levels, the probability that 
this method gives incorrect results (i.e., nv v > to^ or 
to" < to^ 1 ) is also shown in this table, as well as the 
probability of detecting finite m Ve (i.e., m v > 0). These 
results suggest that the estimated confidence levels are 
roughly valid for m^ c = and 3 eV. For m^ LC = 5 eV, 
to^ is systematically smaller than to^ c and 28 trials out 
of 100 MCs give incorrect results of m u < 5 eV. However, 
this systematic error is not greater than 1 eV and we can 
detect finite m Vc with a probability of 99 % if m„ o = 5 
eV. The origin of this systematic error is difficult to un- 
derstand clearly, but probably it is the gradual hardening 
of the neutrino spectrum. The probability of detecting 
finite mass is about 50% for m Ve = 3 eV, and we can 
conclude that this method is marginally sensitive to the 
electron neutrino mass of 3 eV, and can easily detect m Vt , 
of 5 eV. This strategy also gives a fit of neutrino effective 
temperature Tp*, and the average of Tp * is also given in 
the table. This Tp agrees well with the true neutrino 
spectrum, considering that average 9 e energy is 3.15T Pe 
in FD distribution and that of the numerical supernova 
model in this early phase is about 10-12 MeV. (FD-fit 
average energy is a little lower than the true value be- 
cause of the deviation of the true spectrum from the FD 
distribution. See Fig. 9 of ref . ||.) 

Now let us consider the dependence of the proposed 
strategy on the distance to a supernova. There are two 
competing effects: available number of events becomes 
smaller with increasing distance, while the time delay 
due to the finite mass increases. In order to see which 
is more effective, we have tested the proposed method 
against supernovae at D = 5 and 20 kpc, with m^f c = 
eV. Statistical average of to" for 100 MC simulations is 
2.4 ± 0.6 and 3.1 ± 0.7 [eV] for D = 5 and 20 kpc cases, 
respectively (Table 1). Here we have used 400 and 100 
as the values of N cut , which are found to be appropriate 
from TO^-iVcut diagrams. Combined with the fact that 
the average of to" c for D = 10 kpc case is 2.8 ± 0.7 [eV], 



the sensitivity becomes slightly better with decreasing 
distance, but the dependence is very weak and smaller 
than statistical dispersion. Therefore we conclude that 
the sensitivity of the proposed method is roughly the 
same for any collapse-driven supernova in our Galaxy. 
In the above estimate, we have abandoned low-energy 
events below 10 MeV, to avoid 7-ray events induced by 
neutral current reactions of (or v T ) with 16 O which 
are expected to be roughly the same number with v e p 
events in 5-10 MeV If we could somehow remove or 
effectively subtract these noises and apply the method 
to all events with detection energy greater than 5 MeV 
(threshold of the SK), then the average of to" could be 
as low as 1.6 ± 0.4 [eV] for a supernova at D = 10 kpc 
(Table 1). 



IV. DISCUSSION 

Finally we discuss some implications of the reported 
sensitivity of supernova neutrinos to the electron neu- 
trino mass. Currently there are some hints on nonzero 
neutrino masses in astrophysical and cosmological obser- 
vations, and here we consider the following three: the 
solar and atmospheric neutrino anomalies and hot dark 
matter in the universe. The standard cold dark matter 
model of structure formation normalized by the COBE 
data with Qq = 1 is known to be inconsistent with the 
clustering properties of galaxies or clusters of galaxies, 
and the cold+hot dark matter model with f2„ ~ 0.2-0.3 
is one of some possibilities to resolve this discrepancy, 
where VL V is the fraction of hot dark matter in the crit- 
ical density of the universe (e.g., Ref. |]l0|). The most 
promising solution for the solar neutrino problem is the 
MSW solutions with Am 2 ~ 10~ 5 [eV 2 ] JlJ, and the at- 
mospheric neutrino anomaly can be explained by the neu- 
trino oscillation with Am 2 ~ lO^-lO" 1 [eV 2 ] The 
only way to combine these hints without sterile neutri- 
nos is an almost degenerate mass hierarchy of neutrinos 
||, with m v = 4.6(/i/0.7) 2 (fV0.3) [eV] for all the three 
generations of neutrinos, where h is the Hubble constant 
-Ho/(100km/s/Mpc). Oscillations between v e <-» and 

<-> v T give the solutions for the solar and atmospheric 
neutrino problems, respectively. The proposed method to 
constrain TO„ e from supernova neutrino bursts can probe 
m Ue as low as this scale and might detect a finite m Vc if 
the mass hierarchy of neutrinos was actually degenerate. 

Let us briefly discuss effects of possible neutrino os- 
cillations. If the mixing angle between v e and others is 
order unity, the observed v e spectrum would be a mixture 
of original D e and P M (or v T ) due to the vacuum oscilla- 
tion. (The MSW matter oscillation in supernovae is not 
effective for antineutrinos unless the mass hierarchy is in- 
verse.) If the hierarchy of neutrino masses is not degen- 
erate, i.e., toi <C to-2 <C to 3 as observed in charged lep- 
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tons, the mixture of light neutrinos with negligible time 
delay and heavy neutrinos with significant delay would 
make the proposed method inapplicable to measurement 
of m Ft . However, in case of the almost degenerate hi- 
erarchy, the effect of oscillation is only deformation of 
observed P e spectrum due to contamination of original 
C^'s, and the proposed method is still applicable because 
a detailed shape of v e spectrum is not important in this 
method. 
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Figure 1. The result of the proposed strategy of measuring ra v , 
applied to four Monte-Carlo realizations of the SK detection of a 
supernova at the Galactic center with simulated m^f- c of 0, 3, 5, 
and 7 eV. Best fit m u , (filled circles) and 95 % C.L. lower- and 
upper-limits (open circles) are shown as functions of iVcut. 
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TABLE I. The results of application of the proposed method for measuring m v , to virtual data of supernova neutrinos 
produced by 100 Monte-Carlo simulations of the SK detection, where D is the distance to a supernova, N m t the number of 
events used, eth the threshold energy in the analysis, and mil and ml, {ml,) are the best-fit mass and 95 % C.L. lower (upper) 
limits, respectively. 
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